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We investigate chiral and conformal properties of the lattice QCD with eight flavors (Nf — 8) 
through meson spectrum using the Highly Improved Staggered Quark (HISQ) action. We also 
compare our results with those of Nf = 12 and Nf = 4 which we study on the same systematics. 
We find that the decay constant of the pseudoscalar meson "pion" tt is non-zero, with its mass 
Mtt consistent with zero, both in the chiral limit extrapolation of the chiral perturbation theory 
(ChPT). We also measure other quantities which we find are in accord with the tt data results: The 
p meson mass is consistent with non-zero in the chiral limit, and so is the chiral condensate, with 
its value neatly coinciding with that from the Gell-Mann-Oakes-Renner relation in the chiral limit. 
Thus our data for the Nf = 8 QCD are consistent with the spontaneously broken chiral symmetry. 
Remarkably enough, while the Nf — 8 data near the chiral limit are well described by the ChPT, 
those for the relatively large fermion bare mass m/ away from the chiral limit actually exhibit a 
finite-size hyperscaling relation, suggesting a large anomalous dimension 7m ~ 1- This implies that 
there exists a remnant of the infrared conformality, and suggests that a typical technicolor ( "one- 
family model") as modeled by the Nf — 8 QCD can be a walking technicolor theory having an 
approximate scale invariance with large anomalous dimension 7^ ~ 1- 
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I. INTRODUCTION 

The origin of mass is the most urgent issue of the particle physics today. Ahhough the LHC has discovered a 
125 GeV boson roughly consistent with the Standard Model (SM) Higgs boson, there still remain many unsolved 
problems with the SM, which would require physics beyond the SM. One of the candidates for the theory beyond 
the SM towards that problem is the Walking Technicolor (WTC) [l| having a large anomalous dimension 7^ — 1 
and approximate scale invariance due to the almost non-running ("walking") coupling '2*1, which is based on the 
scale-invariant gauge dynamics (ladder Schwinger- Dyson equation [3|,|J]). Actually, WTC predicts a light scalar 
Higgs-like composite, technidilaton, a pseudo Nambu-Goldstone boson of the spontaneously broken approximate scale 
symmetry, which may be identified with the 125 GeV boson [6*]. 

The walking behavior can in fact be realized in the "large Nf QCD" , QCD with large number of (massless) flavors 
Nf, which possesses the Caswell-Banks-Zaks (CBZ) infrared fixed point (IRFP) a* = a^,{Nc^ Nf) (< 00) of the two- 
loop beta function, for iVjr(~ 8) < Nf < ivj^^^ = lliVc/2(= 16.5) in such a way that a* ^- as TV/ ^ ivj^^\ where 

7V}^^) is the maximum number to keep the asymptotic freedom. Due to the CBZ IRFP there exists an approximate 
scale invariance a{^) ~ a* in the infrared region < /i < Aqcd ("infrared conformality" ) , while such a scale symmetry 
is lost for the ultraviolet region /i > Aqcd where the coupling runs as in a usual asymptotically free theory^ When 
Nf is near so that a* is strong enough to trigger the spontaneous chiral symmetry breaking (S^SB), the exact 
IRFP would actually be washed out by the dynamical generation of a quark mass mo 7^ through a continuous 
phase transition ("conformal phase transition" 0), mo = (a, < aa, or Nf > iV" (> Nj)) to niD 7^ (a, > ctcr, 
or (iVjf <)Nf < Nf), in such a way ("Miransky scaling" 13]) that mo ^ Aqcd • exp (-tt/^q;, - a") < Aqcd for 
a* ~ ttcr {Nf ~ Nf'), where acr is the critical coupling for the S^SB and N" the critical number of flavors such that 
a^,{Nc, Nf) = Qfcr- The critical number Nf was estimated as Nf ~ ANc ~ 12 Q by comparing the two-loop value of 
the CBZ IRFP with the critical coupling of the ladder SD equation analysis "s*]: Q!,(iVc, Nf) = acr (= (3C2) = 7r/4). 
Now, for Nf{< Nf) very close to Nf, the dynamical mass mn 0) could be much smaller than the intrinsic scale 
mo ^ Aqcd, in sharp contrast to the usual QCD where niD = C'(Aqcd), so that the approximate conformality 
a{^) ~ a* still remains in the wide infrared region mu < < Aqcd as an impact of the would-be IRFP. Such a 
"remnant of conformality" should appear in low-energy quantities. This is the case for the WTC, with the intrinsic 
scale Aqcd being identified with the "ultraviolet" cutoff A of the WTC usually taken as the Extended Technicolor 
(ETC) scale Aetc, and will be the focus of our interest in this paper. 

Although the above results from the two-loop perturbation combined with the ladder approximation are very 
suggestive, the relevant dynamics is obviously of non-perturbative nature, we would need fully non-perturbative 
studies. Among others the lattice simulations developed in the lattice QCD would be the most powerful tool to 
investigate the walking behavior of the large Nf QCD. Actually, there were some pioneering works on the large Nf 
QCD in somewhat different contexts |10l4l3| . and more recently there have been many lattice studies towards the 
above problemflj]. The immediate issues are: What is the critical number A''"? What is the signatures of the walking 
theory on the lattice? In particular, the above two-loop/ladder studies would suggest that the walking theory if existed 
might be in between Nf = 8 and Nf = 12. As to Nf = 12 there have been many analyses including those of ourselves 
which are consistent with the theory being inside the conformal window [l^, [U, [T5l - [2i | , although some works prefer 
the SxSB phase [2^, [2^. There were also simulations on Nf ~ 10 [l^l consistent with the infrared conformality. We 
thus are interested in Nf = 8 as a candidate for the walking theory. 

Actually, the Nf = 8 is particularly interesting from the model-building point of view [28|: A typical technicolor 
model is the so-called one- family model (Farhi-Susskind model [23] ) which has a one-family of colored and uncolored 
weak-doublets techni-fermions (techni-quarks and techni-leptons) corresponding to each family of the SM quarks and 
leptons. It can embed the technicolor gauge and the gauged three generations of the SM fermions into a single gauge 
group (ETC) and thus is the most straightforward way to accommodate the techni-fermions and the SM fermions into 
a simple scheme to give mass to the SM fermions. Thus if the Nf = 8 turns out to be a walking theory, it would be a 
great message for the phenomenology, which is to be tested by the on-going LHC. Actually, the techni-dilaton [l|, Q 
in the WTC for the one-family model is consistent with the present LHC data for 125 GeV boson in a ladder analysis 
and in holographic estimate @ . 

If AT/ = 8 is a walking theory desired for the WTC, it should be inside the SxSB phase Nf ^ 8 < Nf {mo 7^ 0) 
and at the same time be close to the phase boundary with the conformal window Nf > Nf {mo = 0) such that 
mo Aqcd- Now the lattice simulations we are making contain several scale-symmetry breaking parameters, the 



^ The intrinsic scale Aqcd two-loop level is defined as usual by a renormalization-group-invariant scale parameter Aqcd = M ' 

exp ^— J"^'^^ /^ot) ^^'^^ that = 0, where /3(a) = da/d(\n fi) is the two-loop beta function instead of the one-loop one Q]. 

^ As to the immediate questions about the problem with the S, T parameters, see, for example, discussions in Ref. [soil . 




FIG. 1. Schematic two- loop/ladder picture of the gauge coupUng of the massless large Nf QCD as a walking gauge theory in 
the SxSB phase near the conformal window, mo is the dynamical mass of the fermion generated by the SxSB. The effects of 
the bare mass of the fermion m/ would be qualitatively different depending on the cases: Case 1: nif <^ mo (red dotted line) 
well described by ChPT, and Case 2: mj ^ mo (blue dotted line) well described by the hyper scaling. 



fermion bare mass m/ as well as a finite box and lattice spacing a, which do not exist in the continuum theory we 
are interested in. Among others the fermion bare mass m/ obviously distorts the ideal behavior of the breaking of 
the scale symmetry in a way similar to the continuum theory. Then, disregarding the effects of the lattice parameters 
L and a for the momenl[f|, we may imagine possible effects of the fermion bare mass on the walking coupling of our 
target of study as in Fig.[Tl which is suggested by the two-loop/ladder analysis. 

Case 1. ruf ^ mo <^ Aqcd (red dotted line in Fig. [1]): The chiral perturbation theory should hold in a way similar 
to the real-life QCD with light quarks. 

Case 2. mu <^ mf <C Aqcd (blue dotted line in Fig. [1}: The conformal hyperscaling relation should hold approxi- 
mately with a large anomalous dimension "f„i — 1. 

Actually, the SxSB order parameter to be measured on the lattice is not mo but would be the decay constant of 
the Nambu-Goldstone boson tt extrapolated to the chiral limit: = F^(m/ = 0) which would be expected roughly 
the same as mo'- mo = C){F). 

There is a caveat about the approximate hyperscaling relation to be expected in the Case 2 (mo mf <^ Aqcd )■ 
There are two infrared mass parameters mo and mf which violate the infrared conformality and hence the possible 
hyperscaling relations for the physical mass quantities measured from the spectrum should not be universal but 
do depend on both of them in non-universal ways, in sharp contrast to the hyperscaling relation in the conformal 
window where all the mass parameters from the spectra reflects the deformation by the unique infrared scale- violating 
parameter to/ in a universal way. In particular, when m/ is getting close to the region in Case 1, where tt mass M,r 
and the other quantities such as p mass Mp and -FV behave qualitatively different towards the chiral limit: — ^ 
while the others remain non-zero. 

To date, some groups carried out lattice studies on 8-flavors, with Wilson fermions [I^ililljIHiHl] and with staggered 
fermionM 2^ 5^ l3ll-[33| . The Refs. ilfl, li , 23 , J^] concluded the TV/ = 8 is in the conformal window, but Refs. 
[13, [iH [11, [Hi mi ll^l concluded that the Nf = 8 resides on the chiral broken phase. Even if Nf = 8 is in the chiral 
broken phase, it has not been investigated whether the behavior of this system is QCD like or the walking with the 
large anomalous mass dimension. 

In this paper we study the meson spectrum by simulating the Nf = 8 QCD, based on yet another lattice fermion, 
Highly Improved Staggered Quark (HISQ) [s^, applied to Nf = 8 for the first time. Preliminary reports were given 
in Ref. [36|. HISQ action improves the behavior towards the continuum limit through the improvement of the flavor 
symmetry. The salient feature of our collaboration is that we have been investigating Nf = 4, 8, 12, 16 on the setting 
of HISQ action with the same systematics in order to study the 7V/-dependence of the physics systematically [2l|, [s^ . 
Thus our analyses for Nf — 8 are made in comparison with those for other flavors of our group. 

We first show the data of the meson spectrum, A/^r and F^^, as well as Mp and the chiral condensate {'iptp) for 
/3(= 6/g2) = 3.8 on the x T lattice with and L = 12 - 36 and T = 16 - 48, and m/ = 0.015 - 0.16. We find the 
two regions of to/ having qualitatively different properties: to/ = 0.015 — 0.04 and mf — 0.05 — 0.16. We analyze the 
data based on the Chiral Perturbation Theory (ChPT) [sj, for small to/: m/ = 0.015 — 0.04 (roughly corresponding 



In our simulation we use the parameter region where the effect of the system size is subdominant compared to the mass effect. This 
strategy is different from the one which is advocated by the authors of Refs. [23l. |24|. 
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to Case 1 in Fig. [T] in the above) . We find that the ChPT analysis is self-consistent and find a result consistent with 
non-zero value of F and Mp and vanishing of in the chiral limit extrapolation based on the ChPT (we also estimate 
the effects of the chiral logarithm). The chiral condensate is also non-zero value in the chiral limit extrapolation, 
which neatly coincides with the Gell-Mann-Oakes-Renner (GMOR) relation obtained from the tt data in the chiral 
limit extrapolation. 

As to the large irif {rnj = 0.05 — 0.16) (roughly corresponding to the Case 2 in FigH]), we find the finite-size 
hyperscaling (FSHS) |38l - l4l| holds in this region, when we take into account mass corrections to the FSHS. Note 
that such corrections were sizable [2l| in the large mass region even for Nf = 12, which are consistent to be in the 
conformal window. From the hyperscaling analysis for such a large mj, we find a large anomalous mass-dimension 
7m 1 consistent with that desired by the WTC. This implies that there exists a remnant of the infrared conformality 
where the spontaneous chiral symmetry breaking (SxSB) effects are negligible compared with the mass deformation 
nif. It is the first time that the hyperscaling relation is observed in a theory with S^SB. 

The SxSB feature of Nf = 8 data near the chiral limit are found to be qualitatively similar to those of the Nj — 4 
case: We actually find Nf = 4 data indicate robust signals of SxSB phase. On the other hand, our Nf = 4 data 
indicate no trace of the hyperscaling relation for large m / region in sharp contrast to Nf = 8 data. The Nf = 8 result 
is also contrasted with the Nf — 12 where our previous study concluded that the ChPT analysis with our data was 
not self-consistent, while the FSHS relation held consistently with the infrared conformality. 

This suggests that a typical technicolor ("one-family model") as modeled by the Nf = 8 QCD can be a walking 
technicolor theory having an approximate scale invariance with large anomalous dimension. 

This article is organized as follows: Sec. |ll]presents our lattice simulation setup, calculation of observables, analysis 
method, and the results of the crude analysis of our data. Sec. IIIII shows the analysis based on the ChPT to show 
that Nf = 8 is actually in the SxSB phase. Sec. llVl is to study the remnants of conformality. Sec.|V]is devoted to the 
summary and discussion. Appendices \X\ and IB] summarize detailed numerical results for Nf = 8 and 4, respectively. 
In Appendix [C] we estimate chiral log corrections in iV^ = 8. We analyze FSHS in an alternative method in Appendix 

m. 



II. LATTICE SIMULATION AND THE RESULTS 
A. Lattice setup 

In our simulation, we use the tree-level Symanzik gauge action and the highly improved staggered quark (HISQ) 
action [s^ without the tadpole improvement and the mass correction in the Naik termji^]. It is expected that 
the flavor symmetry in the staggered fermion and the behavior towards the continuum limit are improved by HISQ 
improvement. We carry out the simulation by using the standard Hybrid Monte-Carlo (HMC) algorithm using MILC 
code version 7 14311 with some modifications to suit our needs. One of the modifications is the Hasenbush mass 
preconditioning [44| to reduce the large computational cost of the configuration generation at the smaller m/. We 
measure the mass of the pion M^^, p- meson Mp and the decay constant of the pion F^^ and the chiral condensate (ipip) 
as the basic observable to explore the large- A^/ QCD. 

The simulation in the preliminary report [sg"] , which includes the study of the anomalous dimension, for Nf = 8 is 
carried out at f3{= 6/g^)=3.6, 3.7, 3.8, 3.9, 4.0 for various quark masses and on various lattices, x T, where L is 
the spatial size and T the temporal size. We need to choose as small value of /? as possible to obtain a large enough 
physical volume to minimize the finite-volume effect. From the global survey mentioned above, we found that /? < 3.8 
is too strong to carry out the HMC simulation with HISQ. Therefore, we choose /3 — 3.8 in this article. 

Note that the aspect ratio is kept fixed asT/L = 4/3, in which L — 12, 18, 24, 30 and 36. The boundary condition 
in the spatial direction is the periodic and the one in the temporal direction is the anti-periodic for fermions. We 
take more than 700 trajectories for the ensemble with 4-5 steps for saving the configuration. The error analysis is 
performed with the standard jackknife analysis having the suitable bin size, 40 trajectories. In the following analyses 
the error of the fit result is estimated from the standard deviation of least squares coefficients. See the details of the 
simulation parameter in Tables IXHI IXVH 

We also generate gauge configurations for Nf — 4. From /3 =3.6, 3.7, 3.8 which were investigated in the preliminary 
study [l^, here we focus on /? = 3.7, which is appropriate for our purpose, with high accuracy on 12^ x 18, 16^ x 24 
and 20^ x 30. See Appendix |B] for details. 
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B. Calculation of observables 



We measure the two-point correlation functions of the staggered bihnear pseudoscalar operator which corresponds to 
the Nambu-Goldstone (NG) mode associated with the chiral symmetry of the staggered fermions. The corresponding 
spin- flavor structure is (75 181^5)1 denoted by "PS" in Ref. [45|. The random wall source is used for the quark operator 
for the bilinear, which becomes a noisy estimator of the point bilinear operator with spatial sum at a given time 
slice tQ. We combine quark propagators solved with periodic and antiperiodic boundary conditions in the temporal 
direction (see, e.g., Ref. [i^), which is denoted by "P-I-AP" in this article. In this well-known technique, the temporal 
size is effectively doubled, which enables us to have sufficient range for the fitting. Denoting such a tt correlator as 
Cps{t), this behaves as the following expression in the staggered fermion with P-I-AP prescription: 

Cps W = C (e-^^'* + + Bi-iy , (1) 

where B is the constant term in the oscillation mode and M^r is the mass of NG-pion mode, and C is the amplitude 
relating to the decay. Here we use 

Cps{2t) = Cps(2i)/2 + Cps{2t - l)/4 + Cps(2i + l)/4. (2) 

This linear combination kills the constant oscillation mode, which could originate from the single quark line wrapping 
around the antiperiodic temporal boundary. The mass of NG-pion is obtained by the fit of the two-point correlators 
of Cps from a random source with double period by a fit function with the fit range [tmin,T]^ 

Cps(2i) - C ( e-'"-'' + e-^'^' (21^-2*)) ^ (3) 



where 20 = 0(1 + cosh(M^)). 

The pseudoscalar decay constant, i^^, is obtained through the matrix element of the pseudoscalar operator, 

F. = ^{0\P^iO)K;P), (4) 

by using partially conserved axial current (PCAC) rclatior0. 

We measure Mp from the staggered vector meson operator {'yij4<E)£,i£,4), denoted by PV in Ref. 45|. The asymptotic 
form of the PV correlator at large t may be written as 

Cpv(t) = Ci(e-*^-* + e-A^p^^T-t)) _^ c'2(-l)*(e~*^-* + e-*^--! (^T-t)) 

where corresponds to the mass of the axialvector meson which is the parity partner mode of PV mode in the 
staggered fermion. Since there exists a constant mode due to the wrapping-around effect, we use 

Cpv{2t) = Cpv(2i)/2 -I- Cpv(2i - l)/4 + Cpv(2i + l)/4 . (6) 

Therefore, 

where 2Ci = Ci (1 + cosh(Afp)) and 2C2 = C2 (1 — cosh(A/ai)). Even in the case of Mp ~ Ma^ and Ci ~ C2, we have 
Ci ^ C2 for our typical value of Mp. Then Eq. ([7]) can be approximated to the simple cosh function of the two-point 
correlators of Cpy: 

Cpy{2t) = Ci(e-^^''2* + e-Mp(2T-2t)) ^ (8) 

and we obtain Mp. 

Besides these main channels, we study the masses of mesons interpolated from local operators, a non-NG channel 
(7574 •H' denoted by "SC", and a vector meson (7^ (g) ^i) denoted by "VT" in Ref. [4^, by which we will show 
that the fiavor-symmetry breaking is small in HISQ. These masses are obtained from corner wall source correlator. 



* We use the convention as i<V = v'S/tt, where /,r = 93[MeV] in the real-hfe QCD. 
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FIG. 2. EfTective masses of PS meson, , at L=36. Triangles and other symbols denote results from point sink correlators 
with random wall source and corner wall source, respectively. Fit results with error band obtained from random wall source 
correlator are also plotted by solid lines. 



The effective masses are calculated through C'ps(2t) defined in Eq. Fig. [5] shows typical examples of PS channel 
for the largest volume. The horizontal lines show the results of fitting C'ps(2t) with Eq. ^ with 32 < 2i < 48, where 
plateau is observed. All the fit results are summarized in Tables IXllHXVIl 

We also calculate the chiral condensate {tjjijj) normalized for a single Dirac flavor which can be obtained through 



{■ip{x)ip{x)) = -Tr 



(9) 



where DHisQ{x,y) is the single species (four flavor) staggered Dirac operator for HISQ. Here an average over the 
space-time x is calculated through a stochastic method. 



C. Analysis methods 



We performed the analysis based on the chiral perturbation theory (ChPT) and the (finite-size) hyperscaling, 
as explained in the following; If the system is in the spontaneous chiral symmetry broken (SxSB) phase, physical 
quantities in the spectroscopy, Mh for H = n, p, ■ ■ ■ and F^^, are described by the ChPT. The mass and decay constant 
of TT depend on m/ up to chiral log as 

M^ = C^mf + C^m} + ---, = F + C[mf + C^mj + ■ ■ ■ , (10) 

where F is the value in the chiral limit. 

On the other hand, if the theory is in the conformal window, Mh and F^^ obey the conformal hyperscaling (47| 

Mh oc mj^ , F^ cx mj^ , (11) 

where 7* denotes the mass anomalous dimension 7^ at the infrared fixed point and its value is universal for all 
channels. On the finite volume Mh and F^r are described by the finite size hyperscaling (FSHS) (38l - l4l| on dimension- 
less quantities 

= LMp for p = TT or p, (12) 

or 

^F = LF^, (13) 

given as 

= J-H(Lmf^), (14) 

1 

where H — tt, p or F. The function, J^//, is a some function (unknown a priori) of the scahng variable X = Lmi^'^'' . 
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D. Results 



Spectral quantities, such as Mrr, Mp, F^, are calculated on the gauge field ensembles for the Nf = 8 QCD 

a.t (3 = 3.8, as described in Sect. |lll The m/ dependence of the results is shown in Fig. [3] Large finite size effect is 
observed for smaller nif region on L = 12. 
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FIG. 3. Raw data of observables as a function of m/ for Af^ (top left), (top right), Mp (bottom left) and (tptp) (bottom 
right). 



The expected good flavor symmetry in HISQ action is actually observed in near degeneracy of PS and SC, and of 
PV and VT. See Fig. [H 

Before giving the in-depth analyses in the following sections, let us perform some crude analysis here. Spontaneous 
chiral symmetry breaking leads to non-zero F^r and Mp while vanishing A/^ in the chiral limit. Thus the ratios Ft^/Mt^ 
and Mp/Mj! should diverge in the chiral limit. On the other hand, in the conformal phase the ratios should take a 
constant value near the chiral limit as implied by the hyperscaling relation in Eq. (1111) . 

Now look at Figs.[5]and [Blwhich show that the ratio increases monotonically towards the chiral limit, if one takes 
the largest volume data at each M^. This resembles the Nf = 4 case where S^SB is clearly observed and shows 
clear contrast against the same plot for N j — 12 which are consistent with conformality. This strongly suggests that 
Nf = 8 QCD is in SxSB phase. In order to further study the chiral property of Nf = 8, we carry out ChPT analysis 
in the next section. 



III. CHIRAL PERTURBATION THEORY ANALYSIS 



In order to carry out the ChPT analysis, the finite volume effect has to be taken into account. Fig. [7] shows the 
spatial size L dependence of F^, and Mp plotted from the data on Tables |XlTl IXlIIl IXTVl IXVland IXVll We find 
that the data on the largest two volumes, at least in this ruf range, are consistent with each other. For the lightest 
ruf, since there is only one volume data, we cannot study the finite size effect. We, however, find that the LM^^ in the 
lightest TO/ is bigger than the one of to/ = 0.02 at L = 30 (see. Table IXV)) where the finite size effect is negligible. In 
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FIG. 4. Comparisons of Mt^ and Msc, and of Mp(^pv) and Mp(^vT) as a function of m/ with largest volume data at each m/. 
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FIG. 6. Mp/AU as a function of for Nj = 8 (left) and TV/ = 12 at /3 = 3.7 (right) in Ref. 

the following analysis we understand that there is no finite size effect in the hghtest to/. Therefore, we use the data 
on the largest lattice at each to/ and perform the infinite volume ChPT analysis. 



A. Quadratic fit of Fn 

Let us analyze the behavior of F^, towards the chiral limit. Fig. |5] shows the result of F^^ at each to/. We perform 
the quadratic fit for F^r by varying the fit range of m/. (We will estimate the effect of the chiral log corrections 
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FIG. 7. Ftt (left), Mtt (center) and Alp (right) as functions of L. 



TABLE I. Results of chial fit of with F-,, — F + C'lrrif + C2mj for various fit ranges. 



fit range (my ) 



F 



A'(mf " = 0.015) 



-y('^/ 



xVdof 



dof 



0.015-0.04 
0.015-0.05 
0.015-0.06 
0.015-0.07 
0.015-0.08 
0.015-0.10 
0.015-0.16 



0.0310(13) 

0.0278(8) 

0.0284(6) 

0.0293(5) 

0.0296(4) 

0.0311(3) 

0.0349((2) 



3.74 
4.64 
4.44 
4.18 
4.10 
3.70 
2.94 



11.80 

19.28 

23.2 

26.5 

30.6 

37.0 

54.0 



0.46 
5.56 
4.09 
4.46 
4.06 
7.85 
34.2 



later.) The quadratic fit result of F^^ is written on the TablelH As seen in Fig. [5] and TablelH F is non-zero (~ 0.03). 
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FIG. 8. Results of quadratic fit of for various fit ranges. 

Particularly for the small region, 0.015 < nif < 0.04, the polynomial fit gives the good x^/dof (= 0.46). When we 
include the data at m/ = 0.05, x^/dof jumps up. Although this jump might be caused by the instability due to 
small dof=2, the large x^/dof persists for the range with larger masses, thus, with the value of x^/dof being more 
reliable. This suggests that there is a bound, beyond which the ChPT does not describe the data well, and that 
bound is around mf < 0.05. With this consideration and the good chiral behavior observed for other quantities for 
ruf = 0.015 — 0.04, which we will see below, we chose mj = 0.015 — 0.04 for the fitting range of all quantities. 

For the consistency of the ChPT particularly for the large Nj QCD, the expansion parameter [i^ for the given M^r 
is defined as 

and this quantity is required to satisfy the condition X < 1, which, however, could become easily violated when the 
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TABLE II. Chiral fit of Mp with Mp = + Cfm/ + C^m) for various fit ranges. 



fit range {rrif^ 


Co 


y^/dof 


dof 


0.015-0.04 


0.168(32) 


0.0017 


1 


0.015-0.05 


0.149(33) 


0.098 


2 


0.015-0.06 


0.145(25) 


0.084 


3 


0.015-0.07 


0.144(20) 


0.063 


4 


0.015-0.08 


0.146(16) 


0.052 


5 


0.015-0.10 


0.164(12) 


0.57 


6 


0.015-0.16 


0.189(7) 


1.48 


9 



simulation is made for heavy M^r compared to F. We have X = Oil) in our smallest mj. Thus the ChPT is barely 
self-consistent in contrast to the case of Nf = 12 where A' ~ 40 |2ll|. 

The above analysis suggests that our result in A^j = 8 is consistent with SxSB phase with 

F = 0.0310(13) (16) 
up to chiral log. Effects of the chiral log will be discussed later. 



B. Quadratic fits of Mp and Ml 

Here, we attempt the quadratic fit of A4p and to see whether Mp ^ and = in the chiral limit. 

Fig. [3] and Table HIl are the quadratic fit result of Alp. The chiral limit value of Mp (= Cq) is estimated using the 
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FIG. 9. Results of quadratic fit of Mp for various fit ranges, 
fitting range 0.015 < m/ < 0.04, 

Mp = 0.168(32). (17) 

The left panel on Fig. [TOl shows M^ and the right panel M^/mf as a function of m/. The M^/mf goes to constant 
towards the chiral limit, which is consistent with the leading ChPT behavior. However, the visible slope is observed, 
indicating that there are higher order corrections. This is in contrast to Nf = 4 shown in Fig. [20l We analyze M^ by 
the quadratic fit with the constant term to see whether this constant term becomes zero or not. The result is shown 
in Fig.[TT] In the fitting region 0.015 < m/ < 0.04 the constant term is consistent with zero as presented in Table. [TTTl 
Therefore chiral property of Mp and M^^ is also consistent with that of SxSB. 
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FIG. 11. Quadratic fit of Ml for various fit ranges. 



Chiral condensate 



In this subsection, we analyze the chiral condensate, which is an order parameter of SxSB. We perform a direct 
measurement of the chiral condensate (i/ji/") — Tr [IJ^j^g (x, x)]/4 and compare it with the quantity 



FlMl 

TT TT 



4m/ 



(18) 



which, in the chiral limit, should coincide with the chiral condensate through the Gell-Mann-Oakes-Renner (GMOR) 
relation. Fig. [12] shows the ('0V') ^^id E for each m f . We carry out the quadratic fits for each quantity, whose results 
are summarized in Table IIVI and shown in Fig. 1131 The chiral extrapolations for ('i/'V') and S give good values of 
X^/dof only in the small m/ region 0.015 < m/ < 0.04, though the dof is too small. Both the results in the chiral 



TABLE III. Chiral fit results for Ml with Ml = C J + Clraf + Clm) for various fit ranges. 



fit range (m/) 


CI 


xVdof 


dof 


0.015-0.04 


0.0016(13) 


1.21 


1 


0.015-0.05 


-0.0017(9) 


5.90 


2 


0.015-0.06 


-0.0022(6) 


4.18 


3 


0.015-0.07 


-0.0032(5) 


5.00 


4 


0.015-0.08 


-0.0037(5) 


5.44 


5 


0.015-0.10 


-0.0049(4) 


7.28 


6 


0.015-0.16 


-0.0071(3) 


14.8 


9 
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FIG. 12. (i^il^) (Eq. ((9])) and E (Eq. (I18|) ') as a function of m/. The open symbol represents {^tp) and the filled symbol is E. 
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FIG. 13. Quadratic fits of {ipxp) (left) and E (right). 



limit are non-zero, and are consistent with each other, see Fig. 1141 



0.00059(13). 



(19) 



We also estimate the chiral condensate in the chiral limit by multiplying F in Eq. [16] with the value of Af^ /m / in 
the chiral limit obtained from linear fit in Table IIVI 

F"^ ■ I j = 0.00050(3), (20) 

// m/->-0 



4m 



which is consistent with those from the direct and indirect measurements. 

From the analyses up to chiral log of all the observables, F^, M^r, Mp and the chiral condensate, we find that chiral 
property of A^/ = 8 QCD is consistent with that of SxSB. 



D. Chiral log corrections 



So far we have not included the logarithmic correction in the chiral fits. Here we estimate such effects as systematic 
errors on our previous results. 

The logarithmic m/ dependence is predicted by the next leading order (NLO) ChPT for both the M^/nif and 
F-jr [33, whose formulae are given by 

. , , ^21) 

ruf 



2B {^1 - 

f^^f(i- 



X 

Nfx 



log(x) + C3X 

log(a;) + 042; ) , 



(22) 



13 




FIG. 14. E and {iiip) (left panel) as a function of m/. The small region, < m/ < 0.045 and < S, (t/jt/j) < 0.0006, in the 
left panel is enlarged to the right panel. The quadratic fit curves by using the data in 0.015 < nif < 0.04 are shown. The green 
symbol is the value in Eq. (f20)l . 



TABLE IV. Chiral condensate in the chiral limit: The quadratic fit result of {ipip) and E in various fit ranges, {'tp'tp) = 
C^**^ +C^**^m/ + C^*'''^m|^. E = Co +Cfmf +C2mj. The linear fit of M^/mf = C^*^"'''"^' + Cimy yields the combination 



F^M^/i-imf) T^^Co " /4 in the chiral limit. 



fit range 




















(■mf) 


value 


xVdof 


dof 


value 


xVdof 


dof 


value 


xVdof 


dof 


0.015-0.04 


0.00052(5) 


2.65 


1 


0.00059(13) 


1.11 


1 


2.087(18) 


1.34 


2 


0.00050(3) 


0.015-0.05 


0.00037(3) 


11.5 


2 


0.00015(9) 


11.5 


2 


2.126(14) 


5.22 


3 


0.00041(17) 


0.015-0.06 


0.00037(2) 


7.70 


3 


0.00007(7) 


8.12 


3 


2.151(11) 


6.07 


4 


0.00043(13) 


0.015-0.07 


0.00039(2) 


6.51 


4 


0.00002(6) 


6.56 


4 


2.186(10) 


11.4 


5 


0.00047(11) 


0.015-0.08 


0.00041(2) 


6.23 


5 


-0.00003(6) 


6.22 


5 


2.204(9) 


14.2 


6 


0.00048(9) 


0.015-0.10 


0.00041(2) 


5.19 


6 


-0.00013(4) 


6.75 


6 










0.015-0.16 


0.00056(1) 


18.9 


9 


-0.00026(3) 


7.19 


9 











where the expansion parameter is denoted by a; = 4Bmf /{AttF)'^ , and B,F,C3 and C4 are the low energy constants. 
Our data do not have such logarithmic dependence even in the lightest my region as shown in the previous subsections. 
Actually, such a fit leads to a large x^- This is due to the fact that our m/ is much heavier than the region where 
the NLO ChPT is applicable. The log correction of the F^^, however, is enhanced by the Nf (Eq. (|22p ). so that the F 
might be largely affected by this correction, especially in this large Nf theory. Thus, we attempt to estimate the size 
of the correction by matching our polynomial fit results to the NLO ChPT at m/ such that X = 1, with X defined 
in Eq. where F should read the re-estimated one in this analysis. The details of the analysis are explained in 
Appendix [Cl A reasonable value of the A" < 1 is realized only in the region, m/ < 0.002, much lighter than the m/ 
used in our simulation. From the analysis we find that the log correction reduces the value of F by about 30% from 
the result with quadratic fit. 

The other low energy constants including B are obtained simultaneously. The log correction of the chiral condensate 
is estimated from the GMOR relation, {ipip)\^ — BF^/2, where the values of F and B estimated in this analysis 

are used. We find that the chiral condensate is reduced by roughly half from the result with quadratic fit by the log 
correction. 

Apart from the log correction, we also estimate the systematic error from other sources. While we adopted the 
quadratic chiral fits for F and (V'V') |„^_^q5 hnear fits work with reasonable with the same fitting range. The 
differences are counted as systematic errors. For the chiral condensate, the largest difference from the result with the 
direct measurement to one of the indirect measurements is counted as a systematic error. 

The results for the decay constant and the chiral condensate at the chiral limit in this work are 

^^ = 0.031(l)(i?o), (23) 
W)L^_o = 0-00052(5)(t«g), (24) 

where the first and second errors are statistical and systematic ones, respectively. The lower systematic errors are 
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TABLE V. Power fit results of for various fit ranges, using 
with minimum mass set to the lightest, nif = 0.015, while t 
mf = 0.16. 



fit range (m/) 


Ci 


7 


xVdof 


0.015-0.04 


0.415(7) 


0.988(19) 


14.8 


0.015-0.05 


0.414(5) 


0.991(15) 


9.84 


0.015-0.06 


0.418(4) 


0.979(12) 


7.88 


0.015-0.07 


0.424(3) 


0.963(9) 


7.35 


0.015-0.08 


0.425(3) 


0.961(8) 


6.15 


0.015-0.10 


0.426(2) 


0.958(7) 


5.31 


0.015-0.16 


0.428(1) 


0.952(4) 


3.98 



\ ~ Cimy^^^~'\ The left table shows the results for the ranges 
i right one does those with maximum mass being the heaviest 



fit range (m/) 


Ci 


7 


xVdof 


0.02-0.16 


0.429(1) 


0.947(4) 


2.22 


0.03-0.16 


0.431(1) 


0.942(5) 


1.94 


0.04-0.16 


0.429(2) 


0.950(10) 


1.23 


0.05-0.16 


0.431(2) 


0.941(7) 


0.66 


0.06-0.16 


0.429(2) 


0.948(9) 


0.44 


0.07-0.16 


0.429(3) 


0.950(10) 


0.52 


0.08-0.16 


0.431(3) 


0.939(14) 


0.20 


0.10-0.16 


0.432(4) 


0.934(19) 


0.23 



coming from the log corrections, while the upper ones from the others. 

It would be useful to estimate physical quantities in units of the F, because in the technicolor model the F is related 
to the weak scale, 

y/N^iF/V2 = 246 GeV, (25) 

where Nd is the number of the fcrmion weak doublets as 1 < Nd < Nf/2. From our result, the ratio Mp/F in the 
chiral limit is given as 

=7.7(1.5)(ti^i), (26) 



F/V2 

where the Mp in the chiral limit is the result of the quadratic fit in Eq. (|17p . 

In this analysis we observe the large corrections of the chiral log term in ChPT. In order to reduce the systematic 
error of the chiral extrapolation and to obtain more accurate predictions in this theory, we will need simulations at 
further small mf region on larger volumes. 



IV. STUDY OF REMNANTS OF CONFORMALITY 



In the previous section we showed that the Nf = 8 theory is in the SxSB phase. However, if this theory is near 
the conformal phase boundary, it is expected that some remnants of the conformal symmetry appear in physical 
quantities. 

Here we start with an analysis of -FV from a different point of view. In the conformal phase the F-j^ obeys the 
hyperscaling relation in the infinite volume, Eq. ([TT|) . We perform the power fit F^ = Cimy^^+^^ with various mf 
ranges, where Ci and 7 are free parameters. The numerical results of the power fit are summarized in the Table Ivl 

The power fit does not work in the lightest mf region, 0.015 < mf < 0.04, in which the F^^ is consistent with 
ChPT analysis and the F is non-zero as presented in the previous section. On the other hand, it is remarkable 
that the fit results in the mass range, w/ > 0.05, are consistent with the power behavior, the same way as the 
hyperscaling relation. Furthermore the estimated 7 is stable in the larger mass region (see the right table of Table 
IV)) . the property expected from hyperscaling. This suggests that, although Nf = 8 QCD is in the SxSB phase, there 
exists a remnant of the conformality. Therefore, in this section, we will carry out further in depth analysis, which 
employs the hyperscaling test on the finite volume for F^^ as well as and Mp, to investigate whether the remnant 
of the conformality really persists. 



A. Finite size hyperscaling test 



If the system is in the conformal window, the data on a finite volume is in good agreement with the finite size 
hyperscaling (FSHS) having a universal value of 7 = 7, at IRFP for observables as given in Eq. ([HI) . In general our 
data of Nf = 8 cannot satisfy the FSHS with universal 7 in the whole range oi mf, because we showed that the 
theory is in the SxSB phase as analyzed in Sec. IIIIl However, because of the power behavior in the middle range 
of the fermion mass as mentioned in the above, we carry out the FSHS test in our data to find a remnant of the 
conformality. 
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For this test, we plot the observables, (Eq. ([T3)) ). and (Eq. ([T2|) ). as functions of X = Lmy^^'*'^-' with 
changing the value of 7. Figures [T51 [TBI and [T71 are the results of the FSHS test oi £,f, £,n and for various 7's: The 
data are aligned (collapsing) at around 7 = 1.0, 0.6 and 0.8, respectively. The optimal values of 7 for the observables 
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FIG. 15. plotted as functions of X with 7 = 0.6 (left), 1.0 (center) and 1.4 (right) for the FSHS test. 
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FIG. 16. ,^,r plotted as functions of X with 7 = 0.2 (left), 0.6 (center) and 1.0 (right) for the FSHS test. 
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FIG. 17. plotted as functions of X with 7 = 0.4 (left), 0.8 (center) and 1.2 (right) for the FSHS test. 



are not universal in this estimate, in contrast to Nf = 12, where the alignments was observed with almost universal 
7 [2l|. It is also noted the existence of alignment for each observable is in contrast to Nf = 4, where no alignment is 
observed (see Appendix [B|) . 

Since Figs. [TSl 1161 and 1171 show the good linear behavior, we carry out a linear ht as the leading approximation of 
FSHS, 

= C + C^X, (27) 

for each observable. This formula becomes the hyperscaling, Eq. pip in the infinite volume limit. In Sec. IIIII we saw 
the ChPT fit worked well for the smallest mass region mf < 0.04 for all the observables. On the other hand, we 
already showed the power-like behavior of the _FV for larger masses mj > 0.05. Thus, we restrict ourselves to fit the 
data in m/ > 0.05 in this analysis. To have good linearity we restrict the data in the larger region, > 8. 

Panels in Fig. [TH] are the fit result of FSHS for ^^r, S,f and from the left to the right. The fitting result is given 
in Table IVIl which is consistent with the ^(7) analysis that does not assume the functional form of fitting given in 
the Appendix |DJ In Fig. [TH] data included in the fit are shown as filled symbols, while those with open symbols are 
not included. The figures show the linear fit works well for ^p (x^/dof — 0.66) and £^f (x^/dof — 0.73) , while the fit 
for has xVdof = 2.52. 
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FIG. 18. Linear fits for the FSHS of (left), (center), and Mp (right). The filled symbols are included in the fit, but 
the open symbols are omitted. The fitted region is m/ > 0.05 and > 8. 



TABLE VI. The 7 fitted by the linear ansatz. The fitted region is m/ > 0.05 and > 8. 





7 






xVdof 




0.5668(26) 


0.049(22) 


2.57766(99) 


2.52 




0.9279(79) 


-0.17(10) 


0.4372(38) 


0.73 




0.798(20) 


0.04(19) 


2.779(69) 


0.66 



The larger x^/dof of the fit might be caused by corrections which are not explained by the simple fit form in 
Ea.p7|. To check the existence of the correction, we fit the data of only on two volumes, and slide the range of 
the volumes to investigate the fit range dependence of the 7. The resulting x^/dof, tabulated in Table IVlIl is better 
than the above fit. The results seem to have a tendency to decease the 7 as the volume decreased. The maximum and 
minimum results deviate from each other by more than two standard deviations, and they also differ from the result 
using the four volumes tabulated in Table IVTl This tendency would suggest that there are corrections to the leading 
behavior of the FSHS in Eg. (|27|) . but it is not clear that this tendency comes from only a finite volume effect, because 
the range of the m/ is also changed as the volume. On the other hand, the results for the and in Table IVIll are 
consistent with each result tabulated in Table IVIl so that we do not expect that there are significant corrections in 
these data. 

We will discuss the types of the corrections of Eq. (I?7l) in the next subsection. 

B. FSHS fits with the correction term 

Since Nf — 8 theory is in SxSB phase, FSHS cannot become accurate by approaching to the chiral limit, which is in 
contrast to the Nf — 12 where FSHS does [2l|. Therefore FSHS is only expected for larger mass region, where mass 
corrections may not be negligible. In fact in the last subsection the decreasing tendency of the j{MTr) depending on 
the fit range is seen, which might suggest that there are corrections in the simple FSHS form in Eq. (P7|) . To include 
mass corrections we assume the same fitting forms as in the Nf — 12 case [2l| as, 

= + C^X + Ci'Lmf. (28) 

Since it is hard to determine the exponent a of the correction term when the fit is performed for each observable 
individually, we fix it in our analysis. Among various choices of the a, we take two values: a = 1 and 2. The first 
choice a = 1 is regarded as an mf correction in the heavy region, and the second one a = 2 may be identified as a 
0{a^) discretization effect. 



TABLE VII. The 7 fitted by the linear ansatz using the data on two volumes. 



L 










6 




7 


xVdof 


7 


xVdof 


7 


xVdof 


(30,24) 


0.5864(61) 


1.09 


0.948(18) 


1.19 


0.84(11) 


0.62 


(24,18) 


0.5720(88) 


0.66 


0.934(23) 


0.88 


0.765(58) 


1.25 


(18,12) 


0.5509(54) 


0.49 


0.924(12) 


0.15 


0.809(28) 


1.00 
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TABLE VIII. FSHS fit with a correction term. The fit function: ^ ^ Qf + Cf X + C" Lmf . The fitted region is m/ > 0.05 
and > 8. 



a = 1 


7 






c¥ 


xVdof 


& 


0.791(57) 


-0.004(25) 


1.74(14) 


1.12(20) 


0.66 




0.965(91) 


-0.016(11) 


0.419(44) 


0.026(65) 


0.74 




0.80(25) 


0.003(190) 


2.78(99) 


-0.01(1.30) 


0.73 




a = 2 


7 








xVdof 




0.620(12) 


-0.001(25) 


2.421(35) 


0.98(21) 


0.74 


Cf 


0.941(30) 


-0.016(10) 


0.432(11) 


0.030(72) 


0.74 




0.792(83) 


0.001(190) 


2.80(23) 


-0.1(1.2) 


0.73 



TABLE IX. Simultaneous FSHS fit with a correction term, ^ = C" + C" X + C-fLm'}, where Q is free parameter, but C2 ~ 0. 
The fitted region is m/ > 0.05 and > 8. Degrees of freedom equals to 32. 



7 




a 


xVdof 


0.9292(82) 




0.879(53) 


0.67 








c? 


cf 


^F 
6 


-0.004(25) 
-0.015(10) 
0.01(19) 


1.290(98) 

0.4366(38) 

2.280(64) 


1.538(39) 
0.61(10) 



Using the fit assumptions we fit each observable with the same data region as in the last subsection, ruf > 0.05 and 
Ctt > 8. The results are tabulated in Table IVlIIl The fit results with both a = 1 and 2 of the show the correction 
term actually takes effect (CJ ^ 0), with reasonable x^/dof. Due to the large correction, the 7 of the ^.^^ is largely 
changed from the one without the correction term in Table IVI) especially in the a = 1 case, and the value becomes 
closer to the ones from the other observables. On the other hand, for the and fits, it is found that the correction 
is negligible, and the resulting 7's are consistent with the ones without the correction, presented in Table IVIl as 
expected in the analyses in the last subsection. While in the a = 1 case, we obtain reasonable consistency of the 
7 from the three observables within less than two standard deviations, we cannot exclude the a = 2 fit. Thus, the 
above analyses would suggest 7 = 0.62-0.97 depending on the observables and also the form of the correction term. 

Since we observed that the values of 7 with Eq. for all the observables become closer to each other than those 
without the correction terms, it might be possible to obtain a common value of the 7 from all the observables using 
the fit including the correction. Thus, we perform simultaneous fit using all the observables Mtt, i^Tr, and Mp with a 
common 7. For simplicity, we assume the absence of the statistical correlations between each data of M^r, F^^, and 
Mp. In the fit we do not fix the value of the a, and treat it as a free parameter. It is expected that the corrections are 
small in the and ^p, so that we first carry out a fit omitting the correction term in the ^p. The result is summarized 
in Table HXl This fit works well, and gives a reasonable value of the x^/dof. The resulting a is close to unity. Similar 
value of a is also obtained from a fit without the correction term in the f p as shown in Table [Xj This means that the 
exponent of the correction term is close to unity in our data, while 7's from the two fits are different each other. The 
difference is regarded as the ambiguity in this estimate. It is also possible to carry out a simultaneous fit without the 
corrections in both the and ^p, and fits with the correction terms for all the observables using the fixed a — 1 and 
a = (3 — 27) /(I + 7), because our data prefer a ~ 1 in the above fits. Note that the last one is inspired by the analytic 
expression of the solution of the Schwinger-Dyson equation [4^. Fig. [T9l shows the fit result with a = 1 as a typical 
result of the simultaneous fit. These results are shown in Table IXII and their 7's agree within the above ambiguity. 
Under the assumption that all the observables give a universal 7, we estimate 7 = 0.78-0.93. 

It is noted that a simultaneous fit including the lighter mass with irif > 0.015 in > 6.8 fails with a large 
X^/dof= 3.5 even if the mass correction is included. This is because the chiral property is dictated by SxSB and 
should not be consistent with universal hyperscaling near the chiral limit. 

To summarize, using the fits with the correction term, we estimated the value of 7 from the three observables, and 
obtain 7 = 0.62-0.97 which depends on the observables and the correction term in the fit form. Furthermore we carry 
out simultaneous FSHS fits with the correction term, since a universal 7 would be expected if the theory is very close 
to the conformal phase boundary even in the SxSB phase. The resulting 7 in the simultaneous fits reads 0.78-0.93. 
These estimated 7's would be identified as the mass anomalous dimension in the walking regime. 
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TABLE X. Simultaneous FSHS fit with a correction term, ^ = C,f + C" X + C" Lm"f , where a is free parameter, but = 0. 
The fitted region is nif > 0.05 and > 8. Degrees of freedom equals to 32. 



7 




a 


xVdof 


0.807(18) 




0.949(74) 


0.77 








cf 






-0.001(25) 
-0.011(10) 
0.04(18) 


1.65(11) 

0.516(15) 

2.754(63) 


1.174(66) 
-0.107(19) 



TABLE XL Simultaneous FSHS fit with a correction term, ^ — Cq + Cf^X + C2 Lm" using several choices of a. The fitted 
region is m/ > 0.05 and ^.n- > 8. 



a = 0.889(55) 








Cf 






-0.005(25) 
-0.0275(98) 
0.53(16) 


1.338(96) 

0.4435(36) 

2.476(39) 


1.494(37) 




7 


= 0.9130(76), xVdof = 


= 1.73, dof = 33 






Of = 1 fixed 




c¥ 




Cf 


Cf 
Cp 




-0.014(24) 
-0.012(10) 
0.01(19) 


1.61(10) 

0.484(30) 

2.60(17) 


1.31(15) 
-0.068(44) 
0.25(24) 




7 


= 0.874(25), xVdof = 


0.75, dof = 32 






a = fixed 




c¥ 


Cf 


Cf 


C- 
Cf 

Cp 




0.020(24) 
-0.011(10) 
0.03(19) 


1.52(39) 

0.572(34) 

2.91(30) 


1.17(35) 
-0.158(52 
-0.15(36) 




7 


= 0.775(56), xVdof = 


0.93, dof = 32 





V. SUMMARY AND DISCUSSION 

In search for a candidate for the Walking Technicolor, we have investigated meson spectrum of Nf = 8 QCD by the 
lattice simulations based on the HISQ action for /3 = 6/(7^ = 3.8, and for the fermion bare mass range irif = 0.015—0.16 
depending on the volume size x T with (L, T) = (12, 16), (18, 24), (24, 32), (30, 40), (36, 48). 

We found that the data of F^, Af^ are consistent with the S^SB well described by the ChPT, suggesting that 
F = FTr{mf — 0) 7^ 0, and MT^{mf — 0) = 0, and Mp is also non-vanishing Mp{mf — 0) 7^ in the chiral limit 
extrapolation, the p mass in units of the tt decay constant was determined and shown in Eq. (|26p . We further found 




FIG. 19. Simultaneous FSHS fit in ^^(left), ^F(center) and ^p(right) with a = 1. The filled symbols are included in the fit, but 
the open symbols are omitted. The fitted region is m/ > 0.05 and > 8. The solid curve is the fit result. For a comparison, 
the simultaneous fit result without correction terms is also plotted by the dashed curve, whose /dof = 83. 



19 



that the chiral condensate also has a non-zero value in the chiral limit, which nicely coincides with those from 

the GMOR relation in that limit. In these analysis we used 0.015 < mj < 0.04. 

The salient feature of our collaboration is that we have been investigating Nf — 4, 8, 12, 16 on the setting of HISQ 
action with the same systematics in order to study the A^/- dependence of the physics systematically [2l|, IsB] • Thus 
our analyses for Nf = 8 are made in comparison with those for Nf = 4 and 12 of our group. The qualitative features 
of Nf — 8 near the chiral limit were found to be similar to those of the Nf = 4 case: The Nf = 4 data indicated 
robust signals of SxSB phase. The result was contrasted with that of the Nf = 12 where in previous study the ChPT 
analysis was not self-consistent, while the finite-size hyperscaling (FSHS) relation held consistently with the conformal 
window. 

We then checked whether this SxSB phase is close to the conformal window, having some remnant of the infrared 
conformality. Remarkably enough, in contrast to the data near the chiral limit {nif < 0.04) indicating the S^SB, 
those for the relatively large fermion bare mass ruf > 0.05 away from the chiral limit actually exhibited a FSHS with 
the scaling exponent ~ 0.57 (xVdof = 2.5), -/{F^) ~ 0.93 {x^/dof = 0.7), 7(Mp) ~ 0.80 (xVdof = 0.7). The 

value of 7 is non-universal depending on the observable, with x^/dof for being large compared to the others. The 
existence of FSHS is in contrast to our Nf = 4 data which we showed no trace of the hyperscaling even for large m f 
region and hence no sign of the conformality as in the ordinary QCD. This implies that there exists a remnant of the 
infrared conformality where the SxSB effects are negligible (for a schematic view, see Case 2 in Fig. [1]). 

These hyperscaling relations were obtained for relatively large nif, nif > 0.05, in contrast to the FSHS in conformal 
window, where FSHS becomes arbitrarily accurate in the chiral limit. Therefore, there could exist large mass correc- 
tions on the hyperscaling relation. If we include possible mass corrections on FSHS for each observable, we obtained 
7, 0.62 ^ 7 ^ 0.97. Here the x^/dof for was improved to the level of the others. 

We then performed a simultaneous fit over all the observables based on certain model fitting functions including 
the mass corrections to see if the universality of 7 can be improved by the corrections. To our surprise we found that 
certain model fitting functions in fact yield a universal value of 7, 0.78 ^ 7 0.93, with the variations depending on 
the model fitting function. Since this result is obtained at a single value of /3, it is important to see this feature holds 
in the continuum limit, which will be studied by carrying out simulations with multiple values of /?. Possibility of 
such a large scaling exponent 7 was discussed using Dirac eigen modes in Ref. [33 |. 

The anomalous dimension discussed in the walking technicolor is of course the value in the chiral limit. Lesson 
from the SD equation analysis in the conformal window tells us [i^ that the value of 7 obtained by the hyperscaling 
relation without corrections is an "effective" one which is distorted by the mass corrections. On the other hand, 
the value determined by the fit explicitly incorporating the mass corrections just corresponds to the 7* {jm at the 
infrared fixed point) in the chiral limit and hence is of direct relevance to the walking technicolor. Although the 7™ 
so determined is the value at infrared scale, it coincides with the one discussed in the walking technicolor evaluated 
at the ultraviolet scale, as far as the infrared conformality for the wide scale hierarchy exists. 

Finally, we should comment on the possible light fiavor-singlet scalar meson in Nf = 8 QCD. The walking technicolor 
predicts [3, [1] a light composite Higgs-like scalar boson, the techni-dilaton^s a pseudo NG boson of the approximate 
scale invariance inherent to the walking dynamics. Actually, it was shown [6| that the techni-dilaton is consistent with 
all the current data of the 125 GeV boson discovered at LHC. Then, if the Nf = 8 QCD behaves as a walking theory 
with approximate scale invariance, it would be expected that a light fiavor-singlet scalar composite does exist. These 
studies are currently under way and details will be reported elsewhere. Since the quantum number of such an object 
is the same as that of the scalar glueballs which may also be light, the lattice analyses near infrared conformality 
should be done with great care about the possible mixing with each other. As such we made preliminary studies of 
both fiavor-singlet scalar and scalar glueballs for Nf = 12. We found a hint of a fiavor-singlet scalar bound state 
lighter than tt for Nf = 12 

Summarizing all our analyses we may infer that a typical technicolor ( "one- family model" ) as modeled by the 
Nf = 8 QCD can be a walking technicolor theory having an approximate scale invariance with large anomalous 
dimension 7,„ ^ 1. 
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TABLE XII. Results of the spectra on V = 12^ x 16, with f^in = 10. 



m.f 








Msc 


Mp(pv) 


Mp(vT) 


{4>^) 


0.04 


1224 


0.0622(15) 


0.4181(110) 


0.4397(113) 


0.5574(370) 


0.5389(360) 


0.02167(4) 


0.05 


1284 


0.0735(12) 


0.4844(65) 


0.4987(72) 


0.6108(81) 


0.6157(71) 


0.02704(5) 


0.06 


1224 


0.0904(15) 


0.4681(79) 


0.4664(120) 


0.6372(237) 


0.6343(218) 


0.03269(6) 


0.07 


1264 


0.1030(9) 


0.5091(38) 


0.5168(64) 


0.6910(102) 


0.6882(87) 


0.03802(7) 


0.08 


1264 


0.1144(6) 


0.5352(23) 


0.5439(26) 


0.7093(69) 


0.7031(53) 


0.04328(6) 


0.09 


1300 


0.1222(7) 


0.5686(17) 


0.5774(24) 


0.7478(54) 


0.7449(54) 


0.04826(7) 


0.10 


1300 


0.1302(6) 


0.6033(19) 


0.6116(23) 


0.7886(65) 


0.7866(61) 


0.05319(6) 


0.12 


2500 


0.1442(4) 


0.6694(12) 


0.6760(13) 


0.8556(43) 


0.8517(43) 


0.06265(4) 


U.i4 


ZOUU 


0.1565(3) 


0.7384(11) 


0.7460(13) 


[J.yoZL(oo ) 


U.yoi ( (oy ) 


U.U/ ioH4J 


U.16 


3o24 


0.1676(2) 


0.8056(8) 


0.8142(9) 


l.UU32(29j 


l.UU29(26 j 


U.UoUo9(3j 






TABLE XIII. Results of the spectra on V = 


1 s3 V 94 with f ■ 


— IQ 




nif 




F^ 




Msc 


Mp^pv) 


Mp{VT) 




0.04 


712 


0.0823(5) 


0.3421(29) 


0.3445(33) 


0.4901(96) 


0.4842(107) 


0.02295(3) 


0.05 


752 


0.0910(6) 


0.3886(15) 


0.3908(19) 


0.5323(84) 


0.5248(79) 


0.02818(3) 


0.06 


904 


0.0999(5) 


0.4317(15) 


0.4351(17) 


0.5900(63) 


0.5943(52) 


0.03334(2) 


0.07 


1064 


0.1090(5) 


0.4734(10) 


0.4777(13) 


0.6436(63) 


0.6398(64) 


0.03849(3) 


0.08 


844 


0.1170(5) 


0.5144(10) 


0.5181(12) 


0.6782(49) 


0.6776(56) 


0.04354(4) 


0.10 


876 


0.1315(4) 


0.5948(11) 


0.5993(12) 


0.7729(65) 


0.7698(56) 


0.05334(3) 








Appendix A: 


Data tables in 


Nf = 8 case 







The Nf = % simulations are done at /3 = 3.8 with the fixed aspect ratio L/T = 3/4: (i, T) = (12, 16), (18, 24), (30,_40) 
and (36,48) and with various fermion masses m/. Resultant values of F^r, M-^, Msc, Mlpi^py-^, Mpi^yrp-^ and (V'V') 
for each parameter, together with the number of trajectories A'trj used for the measurements after thermalzation are 
summarized in Tables Kni IXllIl IXTVl |XV]and IXWl 



Appendix B: Simulations of Nf = A QCD 

In this appendix, we show the results of simulations for Nf = 4 QCD. We take (3 = 3.7, and (L, T) = (12, 18), (16, 24) 
and (20,30). For each lattice size, we carry out simulations for various input my, and resultant values of i^^, Mtt, 
Msc and {i'i') for each parameter, together with the number of trajectories A'trj used for the measurements after 
thermalzation, are summarized in Tables IXVIH IXVIIII and IXIXI 

In Fig. [201 plot (left panel) and i^^r as functions of m/. Curves in each figure are obtained by fitting 
ciTUf + C2ni?f and C3 + C/^mf + csm^ to the largest-volume (i.e., L — 20, T — 30) data of and i^^, respectively. 
Fit results are ci = 4.36(3), C2 = 4.3(1.0), C3 = 0.0873(10), C4 = 1.846(87) and cg = -12.3(1.6). These plots show 
typical behavior of a theory which is in the SxSB phase, namely, is well fitted by a linear function (plus small 
quadratic correction) of mj, and i<V clearly has a non-zero value in the chiral limit. As further confirmation, we also 



TABLE XIV. Results of the spectra on V = 24^ x 32 with t^in = 22. 



rrif 


iVtrj 


F^ 


M^ 


Msc 


Mp^pv) 


Mp(VT) 




0.02 


744 


0.0566(8) 


0.2330(25) 


0.2367(37) 


0.3508(117) 


0.3461(111) 


0.01188(3) 


0.03 


728 


0.0715(4) 


0.2832(14) 


0.2851(14) 


0.4044(96) 


0.4018(101) 


0.01746(2) 


0.04 


864 


0.0823(2) 


0.3353(7) 


0.3382(7) 


0.4678(57) 


0.4693(57) 


0.02290(1) 


0.05 


752 


0.0918(5) 


0.3826(10) 


0.3851(11) 


0.5274(54) 


0.5228(53) 


0.02825(2) 


0.06 


1848 


0.1012(3) 


0.4295(6) 


0.4327(6) 


0.5742(58) 


0.5826(55) 


0.03345(2) 


0.07 


864 


0.1088(3) 


0.4731(6) 


0.4767(7) 


0.6288(74) 


0.6345(75) 


0.03853(1) 


0.08 


816 


0.1173(3) 


0.5145(8) 


0.5187(8) 


0.6783(81) 


0.6795(71) 


0.04357(2) 


0.10 


948 


0.1315(3) 


0.5940(5) 


0.5987(6) 


0.7790(65) 


0.7760(68) 


0.05334(1) 
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TABLE XV 


Results of the spectra on V — 


30^ X 40 with 


= 28. 




m.f 






hU 


Msc 









0.02 


996 


n 057S('2"> 


0.2227(9) 


0.2245(10) 


n 3225('75") 


3221 f65~l 


oiisqi"!"! 


0.03 


1044 


0.0709(3) 


0.2801(7) 


0.2818(8) 


0.3967(87) 


0.3940(69) 


0.01744(1) 


0.04 


1060 


0.0826(2) 


0.3354(4) 


0.3377(4) 


0.4718(110) 


0.4730(99) 


0.02294(1) 


0.05 


1108 


0.0918(2) 


0.3834(5) 


0.3859(5) 


0.5317(92) 


0.5302(80) 


0.02822(1) 


0.06 


788 


0.1012(3) 


0.4304(4) 


0.4332(4) 


0.5853(132) 


0.5887(123) 


0.03344(1) 


0.07 


732 


1098^2") 


0.4735(4) 


0.4769(4) 


6349fl41~) 


6329(^120") 


03855 fl~) 








. Results of the spectra onV = 


6b X 4o Wltn tmin 


= 62. 




nif 


iVtrj 




hU 


Msc 








0.015 


1004 


0.0512(3) 


0.1883(5) 


0.1900(7) 


0.2825(107) 


0.2904(71) 


0.00906(1) 


0.02 


808 


0.0579(2) 


0.2216(7) 


0.2229(8) 


0.3202(60) 


0.3232(61) 


0.01189(1) 


0.03 


996 


0.0704(2) 


0.2801(5) 


0.2818(6) 


0.3964(86) 


0.3935(71) 


0.01745(1) 




FIG. 20. (left panel) and Ft^ as functions of m/ for Nf = 4 QCD with /3 — 3.7. Curves in each figure are obtained by 

fitting cim/ + C2mj and C3 + C4?n/ + csm^^ to the largest-volume data of M"^ and F^, respectively. 

plot data of (V'V') obtained from the largest- volume simulation as a function of m/ in Fig. 1211 In the figure, we also 
plotted F^M^/{4mf) (= E) which are calculated from the data of My^ and Fj^ at each m/. Curves in the figure are 
the results of quadratic fits, and resultant values in the chiral limit are 0.00845(14) and 0.00832(21) for and E, 

respectively. 

To estimate the amount of systematic error of chiral extrapolation, we did the same analysis that were done in 
Sec. IIIIDI In the plot, we show the result of ChPT extrapolation which is matched to the quadratic fit result at 
TUf = 0.01, the smallest mass we simulate. The value of F.^- in the chiral limit (F) obtained by this procedure is 
F ~ 0.0730, while the quadratic fit gives F = 0.0873(10). We should note here that there is no visible chiral-log 
behavior in our data in the range of 0.01 < m/ < 0.04, therefore the estimate of the amount of chiral-log effect in 
the chiral limit given here should be understood as the maximum possible. We should also mention here the values 
of chiral expansion parameter X. By using F — 0.0873, the expansion parameter at m"?'" = 0.01 is estimated as 



TABLE XVII. Results of the spectra for iV/ = 4 on 1/ = 12^ x 18. 



nif 


Ntri 


F^ 


M„ 


Msc 




0.01 


500 


0.0858(22) 


0.2373(30) 


0.3022(79) 


0.01471(48) 


0.02 


500 


0.1141(11) 


0.3016(26) 


0.3502(52) 


0.02332(28) 


0.03 


500 


0.12900(69) 


0.3706(11) 


0.4087(28) 


0.03056(15) 


0.04 


500 


0.13842(67) 


0.4283(14) 


0.4705(21) 


0.03683(16) 


0.05 


500 


0.14893(59) 


0.4778(10) 


0.5173(27) 


0.04291(29) 
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TABLE XVIII. Results of the spectra for iV/ = 4 on \/ = 16^ x 24. 



ruf 








Msc 




0.005 


500 


0.08195(76) 


0.16356(87) 


0.2235(65) 


0.01102(21) 


0.01 


500 


0.10258(78) 


0.21390(75) 


0.2763(27) 


0.016413(8) 


0.02 


500 


0.11980(71) 


0.2996(11) 


0.3484(16) 


0.02415(12) 


0.03 


500 


0.13059(59) 


0.36739(75) 


0.4099(20) 


0.03060(11) 


0.04 


500 


0.14124(31) 


0.42524(60) 


0.4659(11) 


0.03697(11) 




TABLE XIX. Results of the spectra for N; 


= 4 on V = 20^ 


X 30. 


ruf 








Msc 


W> 


0.01 


180 


0.10443(40) 


0.20966(68) 


0.2611(29) 


0.01637(5) 


0.02 


380 


0.11955(31) 


0.29873(83) 


0.3413(12) 


0.02396(5) 


0.03 


350 


0.13144(33) 


0.36699(89) 


0.4082(13) 


0.030789(36) 


0.04 


200 


0.14160(24) 


0.42579(86) 


0.4663(13) 


0.037317(39) 
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0.03 

W 

a' 



0.01 



0.01 0.02 0.03 0.04 

FIG. 21. {ijjip} as a function of m/ for A^/ = 4 QCD with /3 = 3.7, L = 20, T = 30. The data points indicated by S are 
calculated through F^M'^/{Amf) with the data of Af^ and F-^ at each mj. Curves in the figure are the results of quadratic fits. 
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FIG. 22. Quadratic fit to the largest volume data of Ft^ and ChPT extrapolation. The value in the chiral limit by ChPT is 
F-n = 0.0730, while the quadratic fit gives F^ = 0.0873(10). 

X ~ 0.3 while the one at m™^^ = 0.04 is A" ~ 1.2. This confirms the consistency of using ChPT analysis and thus, 
we conclude that Nf = 4 QCD is in the SxSB phase. 

In Fig. [221 we plot the values of AfJ^ together with values of A/^ for each m/. This comparison shows the amount 
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of flavor-symmetry-breaking effect in our simulation for Nj = 4 QCD with /3 = 3.7. 
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FIG. 23. Comparison of Mgc and as a function of nif for iV/ = 4 QCD with l3 = 3.7. 

Finally, we show the finite-size hyperscaling test for Nf — 4 QCD by using the data of F-^ obtained here. In Fig.[24l 
we show the finite-size hyperscaling plot for input values of 7 = 0.0, 1.0 and 2.0. As we expect, the data show no 
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FIG. 24. Finite size hyperscaling test of Ftt in Nf = 4 QCD. Input values of 7 are, from left to right panels, 7 — 0.0, 1.0 and 
2.0, respectively. 



alignment in the range of < 7 < 2. This should be regarded as a typical property of QCD-like theory, and contrasted 
to the case of Nf = 8. 



Appendix C: Estimate of chiral log corrections 



In this appendix we estimate the effect of the chiral log correction for the F and (V'V') the chiral limit. Since 
both the NLO ChPT formulae Eqs. ^ and ([22]) contain the B and F, to match the polynomial fit to the NLO 
ChPT formula at the matching point m^, we need to solve the matching conditions using Eqs. (pij) . (|22l) . and their 
derivatives at m^, simultaneously. The m'j dependence for the B and F are plotted in Fig. 1251 For comparison the 
polynomial fit results are shown by the dashed lines. The value of the X in Eg. dTSt evaluated using the obtained F 
at each is presented in Fig. 1261 A reasonable value of A" ^ 1 is obtained only in much smaller m f region than the 
TO/ used in our simulation. At rrij = 0.00199, we obtain X — \ and F — 0.207. 

To study the chiral log correction on the chiral condensate, the value in the chiral limit is also estimated from the 
B and F at each to^ using the GMOR relation, (i/)-!/))!^ = BF^ /2 presented in Fig. [571 When one chooses the 

matching point where X — 1, one obtains roughly a half value of the polynomial fit result of the direct measurement. 

Similar results are obtained from different analyses, such as using the expansion parameter 2M.^/(47rF.n.)^ instead 
of a; in the NLO formulae Eqs. (|2ip and (15^ . or using (tpip) instead of M^/rrif or for matching. Therefore we use 
only the result in the first case to estimate the chiral log corrections in Sec. IIII Dl 
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FIG. 25. Low energy constants B (left) and F (right) estimated by matching NLO ChPT formulae at the matching point mj. 
The dashed and dotted lines denote the polynomial fit result and m'f at A" = 1, respectively. 
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FIG. 26. X estimated by matching NLO ChPT formulae at the matching point ruf. The dotted line denotes X = 1. 




FIG. 27. E|mj->o estimated by matching NLO ChPT formulae at the matching point raj. The dashed and dotted lines denote 
the polynomial fit result and m'^f at X = 1, respectively. 



Appendix D: P{^/) analysis 



In this appendix in order to perform the analysis that does not assume the functional form of fitting in FSHS, we 
consider P(7) defined in Ref. |2]| . 

To quantify the "alignment" we introduce an evaluation function P{j) for an observable p as follows. Suppose 
is a data point of the measured observable p at Xj = Lj ■ mj^^^"''^^ and S^^ is the error of f^. j labels distinction of 
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parameters L and m/. Let be a subset of data points {{xk,^k)} from which we construct a function f''^\x) which 
represents the subset of data. Then, the evaluation function is defined as 



(Dl) 



where L runs through ah the lattice sizes we have, the sum over j is taken for a set of data points which do not belong 
to Kl which includes all the data obtained on the lattice with size L. M denotes the total number of summation. 
Here, we choose for the function Z*-^^' a linear interpolation of the data points of the fixed lattice size L for simplicity, 
which should be a good approximation of ^ for large x. 

This evaluation function takes a smaller value when the data points are more closely collapsed to the line /'■^^■' 
and thus provides a measure of the alignment. ^(7) varies as the choice of parameter 7 and should show a minimum 
at a certain value of 7 when the optimal alignment of data is achieved. We take it as the optimal value of 7. 

We then estimate the uncertainty of the optimal 7 by properly taking account of the statistical fluctuation of 
as well as its effect to the line f^^'^^- For this purpose, we employ the parametric bootstrap method, in which the 
data point is simulated by a random sample generated by Gaussian distribution with the mean ^-^ and the standard 
deviation 5^^ . The distribution of 7 is thus obtained for a large number of these samples, from which the variance of 7 
is estimated. The systematic error associated with the interpolation will be estimated by choosing different functional 
form with linear or quadratic splines as will be discussed subsequently. 

We use the data as the overlapped region sandwiched between m/ = 0.015 on L = 36 and nif = 0.03 on L — 36. 
Fig. [5S]is the result of P{'y) for all M^, F^r and Mp and there are minima of P{'j). The 7 value at the minimum of 
P{7) is written on Table 
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FIG. 28. P{j) in and Mp in Nf = 8 QCD. 



TABLE XX. Optimal 7 to make P{'y) minimum with statistical error. 





A/. 


F^ 




7 


0.593(2) 


0.955(4) 


0.820(20) 
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